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EDITORIAL  NOTE 


' } 

the  attached  report  by  Dr,  Barrar  attempts  to  give  a com- 
plete summary  of  existence  and  uniqueness  theorems  for  the 
Navier-Stokes  equations  for  incompressible  fluids  and  related 

i 

equations.  This  parallels  Dorothy  L B Bernstein's  useful 
summary  of  existence  and  uniqueness  theorems  for  elliptic^ 
hyperbolic,  and  parabolic  equations  in  her  book  Existehqf? 
Theorems  in  Partial  Differencial  Equations,  Only  abotft  ten 
per  cent  of  the  papers  covered  herein  are  listed  in  Misa 
Bernstein's  summary. 

As  only  about  one-tenth  of  the  references  are  listed 
in  Miss  Bernstein's  summary j and  as  no  other  comparable  compil- 
ation seems  available,  I hope  Dr,  Barrar' s report  will  prove 
useful » 

At  the  e&me  time,  as  it  has  not  been  carefully  edited 
or  checked  independently,  persons  doing  original  research  on 
existence  and  uniqueness  theorems  are  advised  to  consult 
original  sources, 


Garrett  Birkhoff 
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I.  NAY  I Eft**  STOKES  EQUATIONS 
A,  Time  Dependent 

The  time  dependent  Navier-Stokes  equations  for  an  incompres« 

\ O % 

sible  viscous  fluid  are 

0u.(x,t)  / aui(x,t)  2 

(Xa)  — — - + uk(x,t)  — — - — - *3/ <7  u^x^t) 


Ox 


x« 


- 3 (x,t)  + gjfxjt) 


[i  * 1 ...  n) 


dx 


(lb) 


0 t ) 


0 . 


d x , 


In  these  equations,  u^(x,  *)  vector  velocity; 

p(x,  t)  the  pressure;  ~g(x,  t)  the  gravitational  force  per  unit 
mass;  and  the  constants  ° and  \J  , the  density  and  kinematic 

-A. 

viscosity  respectively.  The  usual  boundary  conditions  are  u = 0 
on  any  fixed  rigid  surface4, and  in  an  infinite  region  ^i(x,  t) 

— > U at  infinity. 

In  the  special  case  when  l|(x,  t)  is  the  gradient  of  a 
potential  G (e.g.,  when  it  is  the  vector  acceleration  of  gravity) 
by  setting 


(2) 


P = ~ G 


/O 


the  homogeneous  equation 


■^For  a derivation  and  related  references,  see  (25,  Vol.  1,  p.  9'5]. 
Here  and  below,  x = (x.^,  x2,  x^/  will  denote  vector  position, 
and  = (uj,  Ug,  u^),  vector  velocity. 

2x  = ( x j ...  xn) 
x 

Repeated  indices  imply  summation. 

4See  (25,  Vol.  2,  ps  676]. 

^Onlv  the  fixed  boundary  problem  will  be  discussed  in  this  paper. 

For  a.  discussion  of  the  free  boundary  problem,  see  Levy  [42),  Als'q 
only  the  incompressible  case  will  be  treated,  as  no  existence  c:. 
uniqueness  theorems  axe  known  for  the  Navier-Stokes  equations  for 
compressible  fluids,  with  non-zero  viscosity. 


Ufc{x,t) 
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0ux(xst) 

(3a)  — + «k(x,t) 

2 t K 

3^  (x,t) 

(3b)  - 0 

xi 

4$  obtained  * 


2 *■ 


l/<?2u±(x,t)  - (i  = lsean] 
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B.  Two-Dimensional  Case 

In  the  two-dimensional  case  [i.e«,  x * (x^x^)]  if  g = 
curl  {£,  only  7^  * ^ does  not  vanish.  Hence,  by  taking  the 
curl  of  equation  (3a)  yields,  in  this  case, 

(*>  ^(*,0  + ux(x,t)  (x,t )~y\72  (x,t ),  [i  - 1,2] 


whereby  the  pressure  potential  P is  eliminated. 

Also  in  the  two-dimensional  case,  equation  (3b)  permits  the 
introduction  of  a locally  single-valued  stream  function^  y(x,t) 
such  that 


(5) 


±± 

9y 


u2 


gjr 

2 x 


The  use  of  i<  eliminates  the  necessity  for  (>b). 


Then' 


(6) 


and  hence  equation  (4)  may  be  written 


<-  - v2  * , 

8 


9V  v(x.t ) 

w / — 

2 t 


V rv 

2 x 2 y 


. * \ 
v / 


3tKx.  t.) 


9>y  d X 


.a 


y v ■ v^(x,  t ) 


^See  [52,  Chap.  12,  §2,  equation  (1)]. 
'See  (52,  Chap,  12,  §2,  equation  (2)]. 

Q 

See  (52,  Chap,  12,  §2,  equation  (3)V. 


Similar  simplifications  are  possible  in  the  axially  symmetric 


3, 


case,  but  we  shall  not  describe  them. 


C,  Steady  Flow 

A motion  is  called  steady  (permanent,  stationary)  if  all  the 
partial  derivatives  with  respect  to  time  vanish.  In  steady  moticns, 
the  Navier-Stokes  equations  (la)-(lb)  reduce  to 


(8a) 

(8b) 


y^2Uj(x)  » 

'0U1(x) 
'd  Xj 


< ( x ) 
axk 


3 0 • 


+ 


2>P(x) 

3 x^^ 


[i  = 1 ,,,  n] 


The  usual  boundary  conditions  are  that  u7  is  given  (usually 

zero)  as  u (x,  t)  on  a rigid  surface  S,  aod,  in  the  cf.se  of  an 

— ^ — *) 

irifirlto  region,  a approaches  a given  um  at  infinity. 

For  steady  flows,  equation  (7)  reduces  to 


(9)  „<?**(*)  - atisl  ^V(x)  - -l^(x)  . 

9xp  $'x1  Xj  ^Xg 

It  might  be  conjectured  that,  as  t — >>  <»,  the  solution  of 
any  time  dependent  problem  approaches  the  corresponding  stationary 
solution.  Experimental  evidence  supports  this  conjecture  for 
equations  (j>)  and  (8)  if  R is  sufficiently  small^,  but  not10 
for  large  R, 


c 

"R  is  the  Reynolds  number,  defined  as  the  dimensionless  ratio 
R - ~ , where  v is  a typical  velocity  and  L,  a represent?." 
tive  linear  dimension, 

°oce  [5,  pQ  23], 
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II.  PROBLEMS 


A . Existence  and  Uniqueness}  Time-Dependent  Case 

The  mathematical  problems  concerned  with  the  time  dependent 

Navier-Stokes  equation  (l)  are  of  the  initial  value  type  usually 

ll 

associated  with  parabolic  differential  equations  . If  Q is  a 
region  in  En(x^  ...  x^)  space,  the  problem  here  is  to  find  a 
solrtion  u.,(x,  t)  of  equation  (l)  in  Q and  for  0 * t < T, 
s ich  that  uA(x,  O)  = u°(x)  in  Q;  and  uA(x,  t)  » «iS(x.  t), 
where  u°(x)  is  given  in  Q,  and  "u^(x)  is  given  on  S,  when 
x e S,  where  S is  the  boundary  of  Q.  Such  a problem  will  he 
denoted  by  P^.(Q). 

Various  technical  difficulties  arise  depending  on  whether 
Q is  bounded  or  net,  and  on  the  smoothness  of  the  boundary  S 
of  C,  When  Q is  unbounded,  a condition  such  as  uniformity 
at  infinity, 

(a)  lim  u(x,  t)  =*  h°  (t)  , 

| x | — >*> 

or  flnlteness  of  total  energy, 

</*>  IS  uA(x,t)  u^(x,t)dx  - constant 

Q 

' 1 2 n 

must  also  be  included  to  insure  uniqueness  . When  Q = E , the 

initial  values  a^d  condition  (a)  or  ( p ) are  all  that  are  given. 

To  r.ttain  existe.ve  s.nd  uniqueness  theorems,  S and  the  bound cv y 

val‘ c*-:  .'iust  be  dr--- viv  ->d  more  explicitly  than  above.  This  will 

be  eo in  the  tvae’it  of  the  theorems. 


I-"'  ‘ specie!  c-r,  ; of  Boussinesq  flow  [36,  8334a] 
v;  - Of  U = U(z,t;.  ’he  Navier-Stokes  equation  (3) 
*;  :ces  to  the  heat  c-v  ., - on 


^ U 

•o  + 
O ** 


(z,t) 


i . e . , V -• 
actually  re- 


12 


Hence  it  seems  pl^vs- hi  e that  the  orobles;-.-  <--u-  the  time  dopes:  i 

Nav  .er-Stokes  o»...  » • to  of  the  parabolic  type,. 


3 . Ca3e  of  Steady  Flows 

The  mathematical  problems  concerned  with  the  steady  state 

Navier-Stokss  equation  (8)  are  boundary  value  problems  usually  as- 

1*^ 

sociated  with  differential  equations  of  elliptic  type  "(though  see 
(5-4)  below),  If  0 and  S arc  as  above,  the  problem  is  to  find  a 
solution  of  equation  (8)  in  Q,  that  takes  on  preassigned  values 
on  S.  This  will  be  called  problem  P(Q), 


C « Regularity  of  Solutions 

2 

For  Laplace’s  equation  ^7  U = 0 in  n-dimensions,  a solution 

14 

U in  a region  Q is  analytic  In  the  interior  of  Q.  Moreover, 

if  a portion  T of  the  boundary  is  a manifold  of  class  Cn  and  the 

boundary  value  0 assumed  by  U on  T is  of  class  Cn,  then  in 
any  Q-neighborhood  of  T,  U is  of  class  Cn-\  For  the  steady 

state  problems  discussed  in  this  paper,  analogous  theorems  have  been 
proved  in  some  cases  and  are  probably  true  in  all  cases.  Those  that 

have  been  proved  will  be  listed  later  in  this  paper. 

p 

For  the  heat  equation  t?  U(x>  t)  - 3 U(x,  t)/a  t » 0 in  n 

< < 

dimensions,  a solution  in  a region  Q for  t - t - tj^  is  analytic 

in  the  space  variables  and  of  class  C*  in  t for  x in  the  in- 

< Id 

terior  of  Q and  t < t - t^,  Moroever,  if  the  initial  values 
are  of  class.  Cn  in  (s.  Q,  then  the  solution  U is  of  class 


'T 


irrntat ional  steady  flow  in  the  special  case  of  zero  viscosity 


the  equation  for  the  stream  function  actual' y reduces  to  Laplace's 
equation  t7  = 0,  (This  will  be  discussed  further  on  p,  16  . ) hur.-x.- 
it  6 come  plausible  that  v:ho  problems  for  the  steady  state  Navier-Sv  v. 
equation  are  of  the  boundary  value  type, 

l4See  [33,  P.  220]. 

15See  [34], 
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u. 


Cn  in  x and  of  class  Cn^ ' 


x «•  Qj, 


< < 

*»  ' 1 - *1' 


in  t for 

For  the  time-dependent  problems  discussed  in  this  paper  analogous  the« 
orems  have  been  proved  in  some  cases  and  are  probably  true  in  all 
cases.  Those  that  have  been  proved  will  be  listed  later  in  this  paper. 


D . Some  Unexpected  Results 

To  illustrate  some  of  the  pitfalls  that  must  be  avoided  in  deal- 
ing with  these  problems,  the  following  simple  examples  are  given: 

1*  To  obtain  uniqueness  in  unbounded  regions,  the  solution  must 

l8 

remain  bounded.  The  following  example'1  shows  this: 

H(x,  t)  « exp  (-x2/4t) 


is  a non-vanishing  solution  of 

f ^ /N  \ 


Cio) 


3 2 


9* 


19 

a u 
a t 


o , 


for  t > 0 and 


(11)  lim  K(x,  t)  = O,  for  all  x. 

t— >0 

Hence,  this  function  could  be  added  on  tc  any  solution  of  the  linear 
equation  (10)  to  spoil  uniqueness  properties20. 

2,  Physical  intuition  sometimes  fails  in  mathematical  problems. 
Consider  the  equation  for  the  temperature  in  a one-dimensional  in- 
finite rod  with  variable  heat  conductivity 

2 

(12)  -2-H  (x,  t)  - »(x,  t)  -2_£  (x,  t)  + b(x,  t)  (x,  t) 

"dt  c)x‘  £ X 

+ c(x,  t)  U(x,  t)  + b(x,  t) 
a > 0 . 

TJsee  [JJ. 

l8See 

19 

?0See  ^ootri°te  11  > P*  ^ for  the  relevance  of  this  equation. 

Note  that  ^ligt^  H(x,t)  is  not  defined.  Thie  is  discussed  more  com— 

t—  - O 
p3.et.ely  in  [•)]. 


Physically,  it  is  plausible  that  if  the  tempera -are  is  given  through- 
out for  t « t , the  temperature  is  determiiied  throughout  for  t > 

o 

tQ,  This,  however,  is  not  necessarily  true. 

For  equation  (12)  to  have  at  most  one  bounded  solution  in  the 
half-plane  t > 0 approaching  a given  continuous  function  for  t = O, 

(X  * try 

some  condition  such  as  that  ^ a.  ( y , t)*"  'ay  must  diverge  for  - « 

o 21 

must  be  imposed.  The  following  example  shows  this. 

For  the  self-adjoint  equation 

(13)  —2—  u(x,  t)  + ch2x  (ch2xu  ) « O , 

3t  x x 

if  '£  (x)  is  the  Gauss  normal  function 

x 2 

5(x)  * / e"y  /2  dy  , 

—00 

then  for  each  t < T, 


(U)  u(t,  «)-.!+  ± *-)  - + th  *) 

1/2 (T  - t)  /2(t  - tj 

is  a bounded  solution  of  (13)  that  (because  | th  x|  - 1 and 
£(oo)  =3  l)  converges  to  zero  for  all  x as  t — >>  T-, 


3*  The  two-dimensional  steady  state  Navier-Stokes  equation  (9) 

for  the  problem  U = (0,  0)  on  y = 0,  tj  * (0,  0)  on  y = 1 has 

the  unique  solution  U = 0,  However,  with  the  same  boundary  conditions 

the  incompressible  ideal  fluid  equations  have  an  infinite  number  of 

solutions,  (Any  stream  function  iKy)  of  class  C such  that 

■^1  (0)  = -^i(l)  * O and  ^"(y)  = F(^(y))  for  arbitrary  F will  be 
22 

a solution  , ) 


2T 

22 


See  [21,  p.  125], 
See  [52,  Chap.  11, 


Si,  No.  2], 
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III.  CATALOG  OF  THEOREMS 


A Navier-Stokes  Equations 

Leray  [37]  proved  the  existence,  but  not  the  uniqueness  of  sol- 
utions to  problems  P(Q2),  P (P5),  P(E2  - Q2),  and  P(E5  - Q5 ) for 
Navier-Stokes  equation  (8).  ^ He  remarked  that  his  method  does  not 
generalise  to  P(Qn)  and  P(En  - Qa).  2^  More  specifically,  his 

3 

results  in  E are 

Theorem  1.  Let  0^  be  a bounded  region  in  of  „iass 

2R  -$  p 

I + A1  Ka  4”  o € m ! f V» 


-J  P 

Let  a be  of  class  C on  S with 


( — £ _ 

J a n dS  » 0 
S 


Then  there  exists  a regular  solution  u^  of  problem  P(Q^)  for 
equation  (8)  with  u.^  = ai  on  S. 

Theorem  2:  Let  G?  and  cf  be  defined  a3  in  Theorem  1.  Let 

-4  3 

a be  a vector  defined  in  W of  divergence  zero.  Then  there 

exists  a regular  solution  u±  of  the  homogeneous  equation  (8)  with 

u.  = a.  on  S and  lim  u,  = a,.  Moreover, 

1 1 lx! — — 1 1 


, . fJI  & 


> 3V  \ 

£ -Q  3xk 


Pai)(  ^Ui 

& xk  a 


^4 

— ) dXj  dx^,  dx^ 


is  bounded. 


To  prove  theorems  1 and  2,  Leray  used  an  early  form  of  the 


Leray-Schauder  theorem 


Similar  results  hold  in  E . 


In  the  time  dependent  case  for  the  Pt(E  j problem,  Leray 
[37]  proved 


will  denote  a bounded  set  in  En. 

pll  rl  , I p 

technically,  this  is  because  J »)  dV  does  not  exist  for  n > J> , 

^ L a neighborhood  of  every  pointrp  € S be  represented  parametric?  lly 
7- ~ x j ( u , v ) , 1=1, 2, 3,  with  all  D,x.  being  h-Lipschitz  continuous, 
D.  x,  (u,  v,  )-D,x,  (uAvJ  I - k[  tu,  -u^Th  + lv,-v0|h]  and  the 


v -i  neignour nuoa  ox  every  poinx  p € D oe  represeuiea  pararaeu 
t- .. = x , [ u , v ) , 1=1,2, 3,  with  all  D, x.  being  h-Lipschitz  continuous, 
i.o.t  |D1xi(u1y1)-DiXi (u2v2) | - K[tu1-u2]h  + |v1-v2|h]  and  the 

Jacobian  d (x. x^x., ) / (uv ) being  of  rank  two.  Then  Q is  called 
a region  of  class  ^A-h,  If  Q is  a region  of  class  A-h,  and  all 
D exist  and  are  h-Lipschitz  continuous,  then  Q is  said  to  be  of 
c X \ s s R - h4  This  the  notation  used  in  [4  3 ]• 

^For  a statement  and  proof  of  the  theorem,  see  [’ll]. 


9. 


Theorem  3 : Let  u^x^Xg)  and  3 (x^Xg )/d  be  bounded, 

2 2 27  ” 

continuous,  and  belong  to  L (J3  ).  f Let  ? ^(x^,  x^)/.^  xi  - 0. 

2 

Then  there  exists  a unique  solution  u^x^x^t)  of  problem  Pt(E  ) 
for  equation  (3)  for  all  time,  with  u^x^XgO)  = u^x^g)  and 
such  that  u^XjXgt)  and  Su^x^Xgt  )/0Xj  belong  to  L^(E2) 
and  are  bounded. 

The  method  of  proof  is  that  of  successive  approximations. 
Starting  with  an  existence  theorem  for  the  Stokes  equation  (50). 
Leray  constructed  the  sequence 

n 


3u  n(x,t)  2 

(l5a)  i-~ «=  i>V ( x , t ) + 


(15b) 


ft 

n / » 

y u.  tx,t. ) 


ax. 


3 p‘ 

"d  x. 


(x,t)  + (x,t)  [i  = 1,2] 


= 0 , 


n n-l  3U  in’1(x,  t) 

with  gtn  = ukn  1 (x,t)  — 


<?  x. 


n 


He  then  showed  uA  approached  the  desired  solution. 

For  the  Pt(E^)  problem,  Leray  [39]  was  not  able  to  prove  the 

analogue  of  Theorem  3 for  all  time,  but  only  for  a neighborhood  of 

2 2 

t = 0.  Similarly  for  the  P^(Q  ) problem,  when  Q~  is  convex 

and  of  class  B - h,  Leray  [38]  was  able  to  prove  only  a local 

> 

existence  theorem  in  time,  even  for  u^  = 0 on  S for  t - C. 

2 

Dolidze  [10,  14]  obtained  local  existence  theorems  for  Pt(Q  ), 
Pt(Q5),  Pt(52  - Q2)  and  Pt(E5  - Q5);  and  Cseen  [56,  p.  J2]  pre- 

ceded Leray  in  obtaining  a local  existence  theorem  for  problem 

Pt(E^). 

Actually,  Leray  [38,  39]  went  further  than  merely  proving 
local  existence  theorems  by  successive  approximations.  On  the  one 
hand,  he  gave  criteria  under  which  his  local  solutions  could  be  e.v- 


10. 


tended  for  all  time.  When,  for  example  [65]  in  problem  ?t(£' ) the 
given  initial  conditions  ii^x,  0)  have  an  axis  of  symmetry  and  when 


11m  inf 

a 


0 

-o 


0)  ui(x,  0)  dV  = 0 , 


% 


c 


q 


’■'here  q is  the  distance  from  any  point  x to  the  axis  of  symmetry, 

it  follows  from  Leray's  results  that  the  Navier— Stokes  equations 

00  have  a unique  solution  for  all  time. 

On  the  other  hand,  although  Leray  did  net  derive  the  existence 

of  a regular  solution  for  all  time  for  the  problem  P^(E^),  he  did 

derive  the  existenco  (but  not  the  uniqueness)  of  an  irregular  soiu- 

28 

tion  , which  will  now  be  defined. 

Definition : ui(x,  t)  is  an  irregular  solution  of  equation  (J ) 

when  it  satisfies  the  following  three  conditions: 

A)  0 does  not  differ  from  the  time  axis  by  more  than  a set 

of  measure  zero;  where  0 is  described  as  follows.  Define  an 
interval  of  regularity  eeTe  on  the  time  axis,  one  in  the  interior 
of  which  the  vector  ui(x,  t)  is  a regular  solution  of  equation 
(3),  and  one  for  which  this  statement  cannot  be  made  for  any 
interval  containing  ®e^e»  0 is  the  union  of  these  disjoint 

intervals. 

B)  The  function  >1^  u^(x,  t)u1(x,  tjdx^dx^dx^  is  decreas- 

E'  ~ ' - 

ing  on  0 and  the  initial  time  t = 0, 

C)  When  t — > t’,  then  u^(x,  t’)  — > u^(x,  t)  weakly  in 

. 2 

L norm, 

2 

Leray  [38]  obtained  similar  results  for  the  P^(Q  ) problem  in 
a convex  region. 


“uLeray  called  the»  turbulent  solutions,  Puhem  [18,  1st  ser.  Part  2, 
Chap.  3,  ®S]  speaks  of  irregular  solutions.  Oseen  [56,  p,  ?2]  is 
also  concerned  with  them. 


11. 


Leray  [39]  conjectured  that  some  problems  may  have  irregular 


but  not  regular  solutions.  He  suggested  that  if  one  could  find  a 
function  ^ (x)  belonging  to  L (E'')  and  satisfying 

,2...  _/>, ^Vli  3p 


(16a)  - „Clvi  + xk  U_l]  . . £ v,  5Vi 

3X|j  5 x^ 


fa  > 0] 


3xx 


(16b) 


£ v, 


2*i 


= 0 


then  the  initial  values 

(17) 


ui(x,  o)  = 

ysaT 


would  have  no  corresponding  regular  solution  for  all  time,  but 
would  have  the  irregular  solution 

vif2a(T  - t)“1//2x] 


(18a) 

(18b) 


U, 


y2a[T  - t) 


, for  t e [0,  T) 


H _ 0 , for  t - T. 

However,  as  yet  no  one  has  constructed  the  desired  vector 
v1(x). 

Hopf  [30]  developed  a method  for  establishing  "weak"  solutions 
of  equations  (3)®  Although  he  obtained  his  results  much  more 
simply  than  Leray,  and  they  apply  to  any  region  and  any  dimension, 
his  regularity  results  for  "weak"  solutions  are  not  as  strong  as 
Leray's  for  irregular  solutions,  Hopf's  definition  of  a weak  solu- 
tion is  given  below. 


Definition;  Let  G be  any  open  set  in  E 


n-fl 


( • • • x^ , t ) 


-*  2 

space.  Let  N be  the  class  of  vectors  a of  class  C with  com- 
pact  support  in  G.  If  some  vector  u belongs  to  L^(G)  and 


satisfies 


(A) 


JJ 


5 ad 

a t 


dx  dt  + 


SJ* 


dx  dt  + V 


a 


■-*  12. 

((  a an  ; 

JJ  u.dxct  ~ C 

G a*9**/0  1 


-J 


for  all  a belonging  to  N,  with  div  a * Q, 


(B)  1'  i -^JL.  u.  dx  dt  - 0 , 

G Sx1  1 

for  all  h belonging  to  R, 

then  r,  is  called  a weak  solution  of  (3). 

Hopf  showed  that  a weak  solution  of  (3)  in  G which  belonged 
to  C (G)  is  actually  a regular  solution. 

In  connection  with  (A)  and  (B)  above,  it  is  interesting  to 
note  that  Weyl  [7*0  showed  that  if  U belonged  to  L2(G)  and 
UAh  dV  = 0 for  all  h belonging  to  N,  then  u * u*  a.e«, 
where  u*  is  a regular  solution  of  Laplace’s  equation.  In  other 
words,  he  showed  that  c,  vjak  solution  in  this  case  is  actually  a 
regular  solution.  However,  the  corresponding  result  has  not  beer- 
proved  for  equation  (3),  whose  non-linearity  makes  this  a difficult 
problem. 

2 

Another  approach  to  Pt(E  ) for  equation  (3)  is  that  taker  by 
Kampe  de  Feriet  in  [32].  He  showed  that  the  Fourier  transform 


(19) 


(WjW^t)  = — ]Jg  £(x,y,t)  exp[-i(w1x+w2y)]dx  dy 


of  the  vorticity  ^(x,y,t)  associated  with  the  two-dimensional 
flow  of  an  incompressible  fluid  extending  over  the  entire  (x,y)- 
plane,  with  finite  kinetic  energy  and  ^ e L,  transforms  equation 
(4)  into  29 

(20)  x(v1,nz,t)  = -^(Wj  + w2“)z(w1w2t)  + 


JJ, 

E£ 


'®1W2  " 9°w 


2 — ^2~)z  (®i®2^  )z( 


©1  +©2 


^Assuming  there  is  a solution  to  equation  (20)  leads  to  the  problem 
of  determining  a vector  v in^a  given  region  Q,  when  curl  v and 
div  are  given  in  Q,  and  V is  given  on  S.  A proof  of 
uniqueness  can  be  found  in  [52,p.  ^57]  and  a proof  of  existence  in 
f I mi)  | J e 


Bellman  [4],  starting  with  equation  (2C)  and  the  initial  con- 
dition zCwjWg.O)  - 0{wlW2),  derived  by  successive  approximations 

Theorem  4:  If  max  ^(w.w^l  is  sufficiently  small,  there  is 

a solution  of  (20)  which  is  unique  z(wiw2°)  » and  satis- 

30 

fies  the  inequality'' 


< 8 max  |0(wjW2) 

lz(wiV)l  ■ ~ . 'Trrrz 


[1  + ^(wxc  + w2  )t] 


31 

B,  Incompressible,  Non-Viscous  Fluids 

In  the  case  when  R * «,  equation  (1)  reduces  to  the  equations 

32 

of  motion  of  an  ideal  incompressible  fluid 

^u.(x,t)  £?U.(x,t)  i9P(x,  t) 

(2ia)  + u,  (x,t)  “ — — + giVx>^)  (i=l**a 

at  k 3xk 

(x,t) 

(21b)  = * 0 . 

a xi 

r\ 

The  usual  boundary  conditions  are  u * n « 0,  instead  pf  u » O, 
on  any  fixed  rigid  surface. 

In  this  case,  equation  (4)  reduces  to 

■*'  * 3 Ci(x^) 


(22) 


+ ui(x>t) 


3 X, 


* 0 [i  - 1,2]  , 


C>  t C7 

33 

which  means  that  ^ is  constant  along  the  path  of  every  particle  . 

When  R = *>,  equation  (8)  reduces  to 

^ujL(x)  '^P(x) 


(23a) 

(23b) 


uu(x) 


'2  x. 


- + g4 (*)  [i  * 1 .. . n] 


a«i(*) 


o . 


<?  x. 


^Although  this  is  no  improvement  over  Lerayrs  result,  the  approach 
is  interesting. 

•21 

J For  ideal  flow  for  compressible  fluids,  see  [50,  69,  and  8], 

52See  [52,  Chap.  10,  §1,  No.  6], 

5?ibid.,  No.  1], 
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f)  The  divergence  is  zero 

g)  The  normal  component  of  the  velocity  at  the  boundary  is  zero 

with 

h)  The  paths  of  the  moving  point  dx/dt  ■ a,  dy/dt  * v,  / x » a, 
y » b,  at  t ■ 0,  exist}  and  at  each  instant  there  is  a one  to  one, 
continuous  area  preserving  relation  between  (a,  b)  and  (x,  y), 

(a,  b)  e ft,  (x,  y)  e ft 

i)  The  vortex  density  is  constant  along  the  moving  point,  and 
in  each  finite  time  interval  satisfies  a uniform  Lipschitz  condition 
vji.Ka  respect  to  x,  y,  and  t, 

Lichtenstein  [461  earlier  proved  the  following  local  type  of 
•’.'-istence  theorem  for  the  plane, 

p 

Theorem  6;  Let  a solution  of  problem  Pt(E  ) of  the  ideal  ic- 

< < 

expressible  fluid  equations  exist  for  tQ  - t - ti  with  initial 
values  U (x),  Then  if  the  initial  values  UQ(x)  are  changed 

-T  r ^1  n 

slightly^'3  to  U (x),  a unique  solution  to  problem  P (S^)  will 

e':ist  for  the  initial  conditions  IJ  (x)  for  t - t - t, , 

O''  o 1 

Although  this  theorem  does  not  constitute  an  improvement  over 
the  theorems  of  Wolbiner,  Schaeffer  and  Holder,  Maruhn  [49]  extended 
it  to  axially  symmetric  initial  conditions  in  which  is  not  cov- 

ered by  their  two-dimensional  theorems. 

In  the  steady  state  case,  restricted  to  irrotational  motion, 

36 

the  problem  is  reduced  to  the  Neumann  problem  in  potential  theory  , 
The  existence  and  uniqueness  theorems  for  thi3  problem  can  be  found 
in  Kellogg  [33]  and  Gunther  [27], 


55  | U± (x)  - U°(x)|  < e 

|D1  U°(x) 

- D,  Vj(x)  l < € 

M* 

(1 

ro 

56See  [52,  Chap.  11,  S’ 

5,  p.  440], 
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f)  The  divergence  is  zero 

g)  The  normal  component  of  the  velocity  at  the  boundary  is  zero 

with 

h)  The  paths  of  the  moving  point  dx/dt  * u,  dy/dt  * v,/  x «a, 
y « b,  at  t » O,  existj  and  at  each  instant  there  is  a one  to  one, 
continuous  area  preserving  relation  between  (a,  b)  and  (x,  y), 

(a,  b)  c ft,  (x,  y)  e R 

1)  The  vortex  density  is  constant  along  the  moving  point,  and 
in  each  finite  time  interval  satisfies  a uniform  Lipschitz  condition 
with  respect  to  x,  y,  and  t. 

Lichtenstein  [46]  earlier  proved  the  following  local  type  of 

existence  theorem  for  the  plane. 

Theorem  6;  Let  a solution  of  problem  P^(E  ) of  the  ideal  in- 

< < 

compressible  fluid  equations  exist  for  tQ  - t - with  initial 
values  Uo(x),  Then  if  the  initial  values  UQ(x)  are  changed 

-i  j-  p 

slightly^-3  to  U (x),  a unique  solution  to  problem  P. (E  ) will 

° r*.  t 

exist  for  the  initial  conditions  UQ(x)  for  t^  - t - t^; 

AlthoUgh  this  theofdm  doeS  not  constitute  kri  improvement  over 
the  theorems  of  Wolbiner,  Schaeffer  and  Httlder,  Maruhn  [49]  extended 
it  to  axially  symmetric  initial  conditions  in  which  is  not  cov- 

ered by  their  two-dimensional  theorems. 

In  the  steady  state  case,  restricted  to  irrotational  motion, 

56 

the  problem  is  reduced  to  the  Neumann  problem  in  potential  theory  , 
The  existence  and  uniqueness  theorems  for  this  problem  can  be  found 
in  Kellogg  [33]  and  Gunther  [27]. 


35  - U°(x)t  < € |Di  U°(x)  - Dt  U j ( x ) | 

56See  [52,  Chap.  11,  §3,  P.  440], 


< e 


U = i,2] 


16 


G . Burgers1  One -Dimensional  Analogue 

To  gain  better  insight  into  the  Navier-Stokes  equations, 
Burgers  [6  and  7]  and  Hopf  [29]  have  studied  the  one-dimensional 
analogue 

^(xpt)  au(x,,t)  92u(x  ,t) 

(24)  -i—  + u(x.,  t)  1—  - 1/ — V-  • 

* t 

This  is  evidently  a non-linear  parabolic  equation, 

Hopf  [29]  observed  that  the  transformation 

(25)  0 ■ exp  - ^ [ u dx 
or 

(26)  u=-2l/0x/0 


takes  a solution  of  Burgers1  equation  (24)  into  the  heat  equation 


[27) 


- 


and  vice  versa,  (See  also  J.  D,  Cole,  Quart,  Appl,  Hath.  9 (1951) 
225-236.) 

Thus  by  using  known  theorems  for  (27),  Hopf  established 
Theorem  J:  Suppose  that  uQ(x)  is  integrable  in  every 

finite  x-interval  and  that 


(28) 


f s 


(£)d£  = o(x2),  for  j x | large  , 


1 neu 


(29)  u ( X , t ) = 


\ exp  - P(Xj  y,  t)  dy 


i exp  - f(x»  y> 


dv 


where 


(30) 


'(*,  y>  t)  = ■ 


+ 


o 


uo(7?)  d7? 


is  a regular  solution  of  (24)  in  the  half-plane  t > C that 


17. 


satisfies  the  initial  conditions 
ex  . a 

(51)  u(5,  t)  a?  — > \ u0u)  dt,  as  5 , 

for  every  a*  If,  in  addition,  u0(x)  fs  continuous  for  x = a, 
then 

(32)  u(x,  t)  — ■>  uQ(a),  as  x a,  t — > 0. 

A solution  of  (24 ) which  is  regular  in  some  strip  C < t < T 
and  satisfies  (31)  for  each  value  of  the  number  a,  necessarily 
coincides  with  (29)  in  the  strip. 

The  Case  R «o  . If  u.,  (x,  t,  l/)  is  a solution  of  the 
Navier-Stokes  equation  (3)  for  a given  1/  in  a given  region  Q, 
taking  prescribed  values  on  S,  the  boundary  of  Q,  a problem  is 
to  find  what  limit  (if  any)  u^(x,  ) approaches  as  R — > co. 

From  results  obtained  for  equation  (24),  Kopf  [29,  p.  201] 
conjectured  that  in  the  interior  of  Q as  R — > °°,  u-^fx,  t,  J/) 
approaches  a "generalized"  solution  of  the  ideal  incompressible 
equations  (21),^^'^®  The  approach,  however,  cannot  be  uniform  at 
the  boundary  because  a viscous  fluid  adheres  to  the  wall,  and  an 
ideal  fluid  slides  along  it^9, 

D,  The  Boundary  Layer  Equation 

To  study  the  behavior  of  a fluid  in  the  neighborhood  of  the 

boundary  as  R Prandtl  [62]  derived  the  boundary  layer 

ho  a l 
equation  - x 

^This  "generalized"  solution  is  a weak  solution  of  (21)  in  the 
sense  given  on  p.  11  of  this  report, 

^®Wasow  [71]  treats  the  case  lim  \?2u(x,y)  + A = Af(x,y), 

A— >«  9 X 

Oseen  [56,  Part  3]  also  treats  the  linear  case, 

39See  [73.  p.  383]. 

tu  x = F,  y = 577  where  o//^''u  'i/  * 

3 To?’  physical  background,  see  [25,  Vol,  1,  Chap,  2]  and  [17]. 


18, 

(33)  1*  -il±.  - £L±J!*  . fU)  , 

a 7?  <;> V 7ji  an  an  at  ^rr 

when  7}  is  the  straight-line  boundary  of  the  two-dimensional 

region  Q,  and  f{£)  is  given. 

The  boundary  layer  equation  in  the  steady  state  case  may  be 

conveniently  given  in  a form  due  to  von  Mises  [51].  Setting 
2 ,6  2 

z = k - (-=-&)"  as  the  dependent  variable,  and  f and  % as 

-5  ij 

independent  variables,  von  Mises  reduced  (33)  to  the  non-linear 
parabolic  equation 

(34)  — = \/k2-z  4~f  + hU)x  . 

ae,  a* 

This  form  is  very  convenient  for  theoretical  discussion* 

For  the  von  Mises  equation  (34),  Piscounov  [bl]  proved 
Theorem  8:  Let  z(0,  if)  « 0(f)  ^ z(?*  C)  = k(£)  and 

z ( , «)  = O.  Moreover,  let 

a)  0(0)  « k(0),  0(«)  = 0 

b)  k(0)  - 0(f)  for  f > 0 

c)  0(if' ) - k(0)  - cf  lor  some  >■  0,  with  f < 

42 

d)  k(£)  be  continuous  and  non-decreasing  , 

Then  under  these  boundary  conditions,  the  equation 

(jf)  «*,-(*(*)-  *r1/2 

TV  > 

has  a solution.  Moreover,  let 

e)  0"(?)  - 0 for  some  €2  * 0 with  f - £^r 
Then  the  solution  is  unique. 

When  the  boundary  is  an  angle  of  rp  A,  Falkner  and  Skan  [19] 
42 

This  means  a favorable  pressure  gradient,  i.e,,  the  flow  is 
accelerating. 


19. 

reduced  the  boundary  layer  equation  to  an  ordinary  dif f erentiai 
equation.  Weyl  [73]  gave  an  existence  but  not  a uniqueness  theorem 
for  the  differential  equation.  Using  conformal  coordinates  (|,  77) 
with  ds^  = dx^  4-  dx^  = (d£^  + dTj^ ) e(i-}  77),  he  obtained  the 
more  general  form 


(35) 


hU)(k2(0  - (- i)®] 

3 77 


+ -• 


a2* 


97} 


9 

97? 


9 3>  ij/  _ 24 

9T?  3i  dr? 


s 


where  a)  7?  = O is  the  boundary,  b)  h(l*)  = ^ d log  o(%,  o)/dfc, 
c)  3i//3ri  - u,(l*,  7?)  — > k(£)  for  7?  — > ». 


S.  General  Theory  of  Non-Linear  Parabolic  Equations. 

The  only  results  for  non-linear  parabolic  equations  apply  to 
non-singular  equations  for  bounded  regions.  Although  these  results 
cannot  be  directly  applied  to  the  von  Mises  equation,  it  is  con- 
ceivable that  through  some  suitable  limiting  process  they  can. 

Gevrey  [23,  §§28-29],  through  successive  approximations,  ob- 
tained what  are  probably  the  best  results  for  non-linear  parabolic 
equations.  To  indicate  the  scope  of  his  results,  they  are  given 
below  in  a simplified  form.  By  coordinate  transformations  and 
extension  of  boundary  values  to  all  of  (x,t)-space,  the  theorems 
can  be  applied  under  more  general  conditions  than  given  here.. 

Let  T » |x,  y|x  e [0,  1],  y € [O,  1]J  . 

Theorem  9?  Let  f(x,  y,  z,  p)  ^ be  such  that 

(36)  I f (xl^lzlpl ) ” * (X2^2Z2^2  ^ I ” ^ I xi~x2  ^ I y “* y ^ f 

+ |zx  + *2 | + iPx  - P2I 1 

for  (x,  y)  e T,  jz|,  |p|  < N,  Then  there  exists  a unique  solu- 
tion  to  the  problem 

^ p = 3z/3  x,  q = 3z /b  y. 

24.iL 

Albers  [l]  generalized  this  to  n dimensions. 


20. 


-^2  I x = C 

(37)  - — % - — — = f(x,y,z,p)  in  T with  z = 0 on  ] x = 1 , 

3X  3y  {y  = r 

Theorem  10:  Let  g(x,  y,  z,  d,  q)  be  such  that  d"/3y , 


9g/^z,  ^ g/3  p,  g/ci?  q exist  and  satisfy  a condition  similar 

to  (36)  with  the  addition  of  a term  in  |q^  - | , Lot  <z>z/3 q 

> O.  Then  there  exists  a unique  solution  tc  the  differential 
equation 


(38) 


(9"z  S z f*  = r' 

_ = g(x,y,z,p,q)  in  T,  with  z = O om  x = 1 . 

U = r 


D*  & y 

Another  approach  to  the  existence  theorems  for  the  von  Misos 

and  other  non-linear  parabolic  equations  is  through  difference 

equations.  Luckert  [48]  assumed  equation  (351  ) bad  a solution, 

and  showed  how  to  obtain  it  by  difference  equations,  John  [31 ], 

among  others,  proved  the  convergence  of  solutions  of  difference 

4S 

equations  to  solutions  of  parabolic  equations  in  many  cases  . 


F.  Oseen  Equation. 

In  a given  unbounded  region  Q,  equation  (l)  simplifies  to 


46  ^ ui 

the  linear  Oseen  equation  by  replacing  the  term  u,  -=r~— ~ by 

3>  u , k 

1 where  u,_  = lim  ulr(x,t),  "’hit 

. X I >oo 

gives  an  asymptotic  approximation  for  equation  (1)  at  large  dis- 


the  ^.pproxinicit  ion  u«  c-  y *« uva  v u . — a *iu  « 1 > 

k cy)  k 1 1 ^ k 

ki 


tances.  By  rotating  the  x^  axis  parallel  to  the  ufe  vector,  the 
Oseen  equation  reduces  to 


(39) 


^ui(x,t)  _ aui(x,t) 


<9  t 


+ u 


& X, 


yr^u4(x,t)  ) + g (x,t), 

c?  x^ 

[i  = 1 ...  n]; 


as  usual,  we  also  assume  (lb).  In  the  steady  case,  (39)  reduces 
to 


-'See  John’s  [31]  bibliography  on  this  subject 
^6See  [60]. 


PI. 


(40)  >/92ui(x)  + - G1  ""4  ''  ° Bi(x)»  [ 

d X j ^ 1 


i = 1 


n]. 


The  two-dimensional  homogeneous  equation  (4)  for  the  stream 
function  becomes 

(41)  S.  . u,  zLihijL^  , „ vS(x,t)  , 

at  3 Xx 

and  in  the  steady  state  case, 

(42)  wS(x)  = - u _a!iM  . 

a 

Oseen  [57,  58,  59]  treated  the  homogeneous  equations  (39)  and 
(40)  in  E and  EJ  but  always  in  a domain  of  the  form  - C, 
corresponding  to  a flow  extending  to  infinity  on  one  side  of  a 
fixed  streamline.  His  method  is  that  of  integral  equations,  Far^n 
[20],  also  by  integral  equations,  derived 

Theorem  11:  Problems  P(Q^)  and  ) both  have  one 

and  only  one  solution  when  is  of  class  B - h and  when 


ui = ai 

on 

S,  for 

continuous  aA,  and  when 

A. 

For 

P(Q3), 

Je  a • 7)d3  « ? 

B. 

For 

P(E?  - 

Q ?),  U 0 as  1/r, 

The 

same 

result 

2 

applies  in  E , In  both  the  two- and  three 

dimensional  cases  the  solutions  are  of  class  C°°  in  the  interior. 


G,  rokjter-Planck  Equation 

The  Oseen  equation  is  similar  in  form  to  the  Fokker-Planck 


equation 


sra  (x  o 

[ ^]  ? t ai:Tx't; 


with  the  matrix 


positive 


a special  case  of  (43),  if  V w = 


variable.  Equation  (43)  applies  to 


3 

+ Z b.(x,t)  -j-—  + c(x,t)u+f(x,t ) 

1 <_V  a , 

1 

definite.  Thus  (41)  is  seen  as 
Is  taken  as  the  dependent 
diffusion  with  variable,  non- 


22. 

isotropic  dif fusivity . 

For  the  Fokker-Planck  equation,  Feller  [21]  defined  a funda- 
mental solution  u(x,  t;  J.  ) for  the  homogeneous  equation 


(12')  ut(x,t)  = a(x,t)uxJ£  + b(x,t)ux(x,t)  + c(x,t)u, 


as  follows; 

1.  As  a function  of  x, 

2.  As  a function  of  f, 

adjoint  equation  u t (x,t)  = 

b 

3.  lim  [ u(x,t;^,r) 


t;  u(x,t;  (.')  satisfies  (12*  ) 

L ; u(x,t;  C ) satisfies  the 
(a(x,t)u  ) - (bu)  + cu 

A A A 

I f { x ) a < x < b 
f(i)  d?  - 1 

0 x ■<  a or  b < x 


for  f(x)  continuous,  (if  tho  range  of  x is  not  finite,  some 
additional  condition  on  f(x)  is  needed.  A sufficient  condition 
is  that  f(x)  is  bounded.) 

4,  The  above  properties  hold  throughout  the  infinite  range 
- eo  < x < oo,  tQ  < t < t - tx  , 

and 

CO 

c 

j u(x,t,s,  )u(s,t;£,  L " )ds  = u(x,t;  C)  . 

—00 

Under  the  assumptions 

A’.  For  tQ  3 t - tl  and  all  x;  a,  at,  ax,  axj£,  b,  bx,  c 

are  7-Lipschitz  continuous  in  x and  t. 

B1.  a,  l/a , A , A„,  c are  bounded  whore  A = b - a /2  + yr\  0 
rx  X t 

and  where  0(x,t)  =\  a(y,t)"’1///2  dy, 

' o 


he  proved  in  [21] 
Dressol  [15,  16] 
ently  Weber  [72] 


V 

l~  a 
ij 


ii 
A *#• 


a 

9 


tho  existence  of  a fundamental  solution  for  (12*). 
generalized  this  result  to  equation  (43).  Rec~ 


extended 


v v 

"*i  J 


it  to  the  equation 
n 


/ ( V 

/ — \ A.f 

i- 1 3 


+ aU  + -_\.V  - o. 


2:5. 


Fundamental  solutions  are  used  to  prove  the  existence  of  sol- 
utions of  initial  value  problems,  Using  the  concept  of  a funda- 
mental solution,  Feller  [21]  showed  that  Pt(E^)  for  equation  (12) 
has  a unique  bounded  solution  U that  approaches  a continuous 
0(x)  for  t ~ > Q,  when  f(x,  t)  has  continuous  bounded  first 

47 

derivatives,  and  the  coefficients  satisfy  assumptions  A',  B'  above.1 

When  the  coefficients  satisfy  7-Lipschitz  conditions,  problem 
Pt(Qn)  for  equation  (43)  has  one  and  only  one  solution  for  contin- 
uous initial  and  boundary  values  as  shown  in  Barrar  [3]. 

In  order  to  treat  singular  Fokker-Planck  equations,  for  ex- 
ample, when  a(r^)  = 0 in 

(44)  ut(t,x)  « a(x)uxx(t,x)  + b(x)ux(t,  x) 
and  its  adjoint 

(^5)  vt(t,x)  * (a(x ) v(t,x))  - b(x ) v(t,  x)] 

in  the  interval  [ria  r2].  Feller  [22]  used  an  approach  which 

4S 

depends  on  the  Hille-Yosida  theorem  for  generators  of  semi-groups 


hY 

’’Recently  John  [31*  P»  X 65 ] extended  this  theorem  to  apply  when 
0(x ) is  bounded  and  Riemann  integrable, 

48 

Hille-Yosida  theorem:  Let  S be  a contraction  semi-group  with  range 

dense  in  the  Banach  space  X in  which  positive  elements  are  defined, 
[i.e,,  S is  a set  of  transformations  Tt  (with  t > 0)  such  that  to 
each  x e X and  each  t > C,  there  corresponds  an  element  T.x  e X 
such  that  Tt+gx  = Tt(Tgx)»  Furthermore,  x 

a'  I lTt+h  X " TtXl I ~">  0 as  h — > 0 
*>.  I lTtx!  I - I M I 

c,  x - 0 implies  Ttx  - C] 

The  infinitesimal  generator  A is  an  additive  operator  whose  domain 
is  dense  in  X and  such  that  to  each  x e X and  each  A > 0 there 
exists  a unique  y^  belonging  to  domain  of  A with 

A,  Ay^  - AyA  = x 

B,  A M^i  5 II  x || 

C,  - 0 whenever  x - 0 

Conversely,  an  additive  operator  with  these  properties  is  the  in- 
finitesimal generator  of  a contraction  semi-group  with  range  dense 
in  X,  The  inverse  ye  = (A  I - A)’*x  is  called  the  resolvent  of  S, 

For  a proof,  see  [76|* 


assumed  a'(x).  b(x)  are  continuous,  but  not  necessarily 


24. 


bounded,  in  the  open  interval  (r\v2^  and  a * 0.  To  express  his 
results,  he  defined 

w(x)  = exp  [- 

v/here  xo  € (r,rg)  is  fixed,.  Fe  used  the  regularity  condition  that 
the  boundary  is  called  regular  if  w(x)  belongs  to  the  space 

of  integrable  functions  on  (x  ,r  . )o£(x  ,r  . ) , and  a“i (x)w~* (x ) 

o j o J 

belongs  to  ‘^C(xo,rj),  Thus,  for  example,  for  the  equation  ut  = 
uxx,  Tj  Is  regular  if  and  only  if  r^  is  finite.  With  this  nota- 
tion, he  derived 

Theorem  12:  Under  the  above  conditions,  when  none  of  the 

boundaries  is  regular,  there  exists  one  and  only  one  fundamental 
solution  common  to  (44)  and  (45),.  When  one  or  both  boundaries  are 
regular,  then  a necessary  and  sufficient  condition  that  (44)  and 

(45)  have  a common  fundamental  solution  is  that  there  exist  con* 
stants  p^,  qj  such  that 


x 

^ b(s)  a‘1(s)  da]. 


(46)  q,  lira  u(t,x)  + p.(-l)^  lim  w 1(x)u  (t,x)  - 0 

J x- -> r , ^ x--*r  . 

J J 

(47)  q lim  w(x)a(x)v(t,x)  + p.(-l)^  lim  ( [a(x)v(t,x) ] 

J x-->r j J x — >rj 

- b(x)  v(t,x))  = 0, 

respectively  at  the  regular  boundaries. 

Feller  called  (46)  and  (47)  the  "generalized  classical  boundary 
conditions,”  Pis  method  of  proof  works  only  when  the  coefficients 
in  (44)  and  (45)  do  not  depend  on  the  time, 

Gevrey  [24]  obtained  existence  theorems  fer  P^(QU)  for  the 


parabolic  system  of  equations 

(*8) 


n 

V" 

/ 

k S>su,. 

+ b 

\ 

+ xl- 

j 

fj  A-J  C?  A j 

+ Dv 

i ^ *_ 

c?  t 

hi 

^ U.  V 

-3*7  + \-el\  " f* 


who; 


k 


is  of  class 


l.  , is  positive  definite,  b.  0 and  when 
ij  *- 

B - h and  the  coefficients  are  y-Lipschitz  continuous.  However, 
it  does  not  seem  possible  to  apply  these  restilts  to  Osoen's 
equation. 

At  this  point  it  is  interesting  to  note  the  following 

Theorem  13 s Let  (a^j)  ^or  k = 1 ...  n of  (48)  be  a sym- 

n ]' 

metric  positive  definite  form  in  Q for  0 - t 3 t»  Let  b < 0. 

If  V b 

Let  bhi,  a-M*  ch  9 fu  °f  (48)  all  bo  Y-Lipschitz  continuous  in 
n < •< 

Q for  0 - t - T,  Let  Q be  of  class  B - h.  Then  there  exists 

one  and  only  one  regular  solution  to  problem  ^t^n^  *or  ^ - t 

< 4o 

- T,  for  preassigned  values  of  U^j  if  and  only  if  the  homogeneous 
problem^  has  only  the  trivial  solution  = 0,  k = 1 «•*&  . 

Proof  r^-  The  theorem  will  be  proved,  without  loss  of  generality 

for  Uk(x,t.)  = 0 for  x e S and  t = C.  Let  H be  the  Banach 

-T 

space  of  Vvectors  = (f,  . , . y ) such  that  each  i/.  is  y-Lip- 

^ & Q.. 


schitz  continuous,  with 


\i' 


n 


!^il!T'i*  52  ?or  any  given 


a — r — r 

f e H,  let  U(^)  = (U^  Un)  be  the  vector  such  that 


. k 

aiJ 


3 


2„k 


•3  IT 


k 


ij  J Xj  3 t 


= UJ'(x,t ) = 


x e S 
t = C 


[ k » 1 ...  n ] 


Then  let 


and 


L(U{/)), 


/ 


, k ^ , > k I? 

- . bhi  T + t—  h Uh 

ni  **  x 

' 1 


— ^ 

F=  (f-j  ...  f n ) 
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i.e.,  Uk(x,0)  = U°(x);  Uk(*,t)  - U^Xjt)  for  x £ S,  with 
U°(x)  and  U1(x,t)  given  continuous  functions  thatagree  when 
t = 0. 

i.e.,  U?(x)  = U*(x,t)  = 0 in  footnote  49  above,  and  f,  ^ 0, 


k = 1 ...  n , 

'^The  proof  is  in  the  spirit  of  Schauder  [69,  Chap.  41, 

J^\  I l^11  = *ax  |tf.  | + C„  , where  C.v  = g.l.b.  of  call  Crs  such  that 

T r\  1 1 V 


Q 


n 


:■■!/.  (x,  ,t,  )-tf,  (x0,t0)  j - C[  Ix.-x^l^+l  ''  i and  wh or 

111  -I  (Z.  1 1 d.  I 


A ^ ^ t V* 


Q 


Using  this  notation,  equation  (48)  may  be  written 

V'  + L(U(i/3  = F,  yli,  F e H . 

From  results  in  [>],  it  follows  that  L(U(y))  is  a completely 

continuous  linear  operator  taking  H into  H.  Thus  by  the 

Hildebrandt-Riesz  theorem*^,  equation  (48)  has  a solution  if  and 

S4 

only  if  the  homogeneous  equation  has  only  the  trivial  solution  . 


H.  Stokes's  Equations  (Case  R 4/  C) 

Equation  (l)  may  be  plausibly  simplified  at  low  Reynolds 
number  R « 1 by  assuming  that  all  velocities  are  very  small , 
and  thus  completely  neglecting  the  term  2>u^,/3Xj?  quadratic 

in  the  velocities*  The  resulting  Stokes  equation  ^ 


(49a; 

(49b) 


3 u . (x, t ) 2 p 

) = l/v>  u (x,t)  + — £ (x,t)  + gi(x,t) 

3 t 3 x . 


? t 

3u1(x.t ) 

*5^7 


fi  = 1 .. 


= o 


or,  in  the  steady  state, 


(50a)  j 2uA(x)  + 5 p/<3  xi=Pi(x>t)  [ i “ 1 • • • n] 
(50b)  3ui(x)/3xi  = 0 


is  mathematically  a special  case  of  the  Osoen  equation. 

The  two-dimensional  homogeneous  equation  for  the  stream 
function  becomes 


(51) 


^ *7  *~'{r  ( x , t ) _74  ( A x 

ra — i L - yxy  !>(X,t)  , 


<r? t 


In  the  steady  state,  it  becomes  the  biharmonic  equation 

(52)  = 0 

— 

For  a statement  and  proof,  see  [2,  n,  150], 
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No  known  uniqueness  theorem  is  applicable  to  equation  (43). 
one  in  [241  uses  integral  equations, 

^ _ r t-O  i 

-^oee  [ > >u  j . 


The 


27 


For  the  Stokes  equations,  Odqvist  [54]  derived  a result  sim- 
ilar to  Theorem  11  for  P(Q^),  P(S5  - Vp ) , p(Q2)  and  P(E2  - Q2) 
under  the  somewhat  weaker  restriction  that  Q is  of  class  Ah. 

His  results  also  cover  the  inhomogeneous  equation  when  g^  is 
bounded  and  “y-Lipschitz  continuous,  Odqvist  was  able  to  char- 
acterize the  behavior  of  the  solutions  of  the  Stokes  equation  more 
precisely  than  Faxen  could  characterize  the  behavior  of  the  sol- 
utions of  the  Oseen  equation.  For  example,  Odqvist  proved  that  if 

2 

Q is  of  class  Bh,  and  the  boundary  values  are  of  class  C , 

2 

then  the  solution  is  of  class  C in  Q, 

Leray  [27*  38]  discussed  Pt(E2)  and  Pt(Ep)„  For  example, 
in  [39]  he  derived 

Q ~Z 

Theorem  14  i Let  u^(x)  be  continuous  and  belong  to  L^E'j. 
Let  there  exist  a continuous  function  f(t)  for  t e [O,  t]  such 
that  ( { ; 

t)g1(x,t)dv  - f(t)  . 

E 


Let  3 gi(x,t)/^> exist.  Then  Pt(E^)  for  Stokes  equation  (49) 
has  one  and  only  one  solution  ui(x,t)  with  ui(x,0)  = Uj(x),  and 
such  that  there  exists  a continuous  function  g(t)  defined  for 
t e [0,  T]  such  that 

\\'^  ui(x,t)ui(x,t ) - g(t)  . 


\ 

0 


O 

He  derived  a similar  result  in  [38]  for  Pt(Q^)  with  Q 

convex  and  of  class  B - h,  and  with  u^  = O on  3,  In  [40]  he 

2 

gave  a method  (without  details)  to  transfer  P^(Q  ) and 

x v 

O 2 

?t(E  - Q ) to  Voltorra  integral  equations,  and  thus  gave  a way 
to  obtain  existence  theorems. 

For  pt(7)  and  Pt(E^  - Cp),  Odqvist  [55],  by  a Laplace 
trinnsf orraation,  orthogonal  series,  and  integral  equations,  ob- 

tainod  The  corresponding  existence  theorems  when  Q" 


of  class 


28, 

h 

B , However,  Odqvist,  due  to  technical  difficulties,  had  to  assume 
Dk^ui(x’  °))“  0 , k - O,  1,  *.♦  4 for  x e S, 

Dolidze  [11,  12,  13 ] was  also  concerned  with  Pt  (<P), 

P^ ( fP  ~ Cp  ) , P^(Q^)  and  Pt(E^  - °^ ) • The  existence  theorems 

in  [13]  may  possibly  be  more  general  than  those  in  Leray  [38]  and 
Odqvist  [55]*  However,  J.  Kravtchenko  [Math.  Rev.  9 (194£)  116-117 ] 
has  said  that  Dolidze's  reasoning  is  very  condensed  and  his  hypo- 
theses are  not  clear.  Kneale  [35 1 treated  the  axially  symmetric 
case,  deducing  his  existence  theorems  from  known  ones  for  the  bi- 
harmonic  equation. 

Uniqueness  theorems  for  the  Sr.okes  equations  are  readily  proved. 
The  proof  for  the  steady  state  case  is  given  in  Lichtenstein  [47, 
p.  594];  the  proof  for  the  time  dependent  case  , in  Leray  [39,  P. 

215]55, 

For  the  biharmonic  equation,  Schrftder  [66]  obtained  existence 
and  uniqueness  theorems  for  = 0 for  P(Q^),  P (Cp).  P(E‘'-Q^), 

and  P(E^  - (p ) when  Q is  of  class  Ah,  with  df/,y 

*2  T 

on  S,  For  the  P(Br  - Or  ) existence  and  uniqueness  theorem,  he 

2 2 

postulated  that  £ "ty/  0 vanished  as  1/r.  For  the  P(E  - Q ) 

existence  and  uniqueness  theorem,  he  postulated  that  c?  $/o  was 
bounded  at  infinity.  This  is  connected  with  the  Stokes  paradox, 
which  can  be  stated  as  follows, 

^Sinee  the  Navier-Stokes  equations  are  not  linear,  a more  complicated 
proof  is  needed  for  them.  However,  it  is  very  similar  to  the  proof 
above.  See  [39,  p.  221],  Leray1 s proof  also  works  for  Oseen’s 
equation. 

-^Although  the  proof  on  p,  25  docs  not  apply  to  the  time  dependent 
Stokes  aid  Oseen  equations,  it  may  be  conjectured  that  it  can  be 
modified  to  do  so.  Since  uniqueness  theorems  are  known  for  the 
Stokes  and  Oseen  equations,  the  proof  would  then  yield  an  exist- 
ence theorem.  It  would  also  yield  stronger  bounds  on  the  behavior 
of  the  solution  than  previously  proved,  because  with  this  proof  it 
is  possible  to  go  directly  from  the  Fokkor-Planck  to  the  Cseon  and 
Stokes  equations.  There  is  a similar  relationship  between  elliptic 
equations  and  the  steady  state  Stokes  and  Oseen  equations. 


Although  a solution  of  the  steady  state  homogeneous  Stokes 
equation  (50)  exists  for  the  boundary  conditions  U = 0 on  the 
surface  of  a sphere  and  U — > at  infinity,  the  correspond  nig 

solution  for  the  two-dimensional  problem  of  U **  0 on  the  surface 
of  an  infinite  circular  cylinder  does  not  exist^. 

Odqvist  [54]  also  obtained  existence  theorems  for  the  two- 
dimensional  biharmonic  equation,  but  his  results  are  not  as  gen- 
eral, nor  does  he  describe  the  behavior  at  infinity  as  specifically 
as  Schrfider. 

In  Schr6der  [67],  properties  of  solutions  of  \? ^ = O are 
studied  very  carefully*  For  example,  he  proves  that  if  S is  a 
manifold  of  class  Cn,  and  the  boundary  values  are  of  class  Cn, 
then  $ is  of  class  Cn  in  Other  properties  for  them  can  be 

found  in  Nieoleseo  [55]*  In  [53,  §8],  for  example,  he  gives  the 
proof  that  a biharmonic  (in  fact,  a polyharmonic)  function  defined 
in  Q is  analytic  in  the  intorior  of  Q. 

p 

From  known  results  about  Poisson's  equation  v7  U = g,  it 
is  possible  to  construct  a particular  solution  of  \/rU  = f,-° 

Thus  because  of  t}»®  linearity  of  the  biharmonic  equation,  all  ex- 
istence theorems  for  the  homogeneous  equation  also  apply  to  the  in- 
homogeneous. 

Similarly,  just  as  particular  solutions  of  Poisson’s  equation 
in  infinite  space  may  be  constructed  with  l/r,  particular  solu- 
tions of  Oscen's  and  Stokes s equations  may  be  constructed  with  the 
tensors  given  in  [56],  so  that  the  results  for  the  homogeneous 
equation  again  apply  to  the  inhomogeneous. 


^^See 

^SLet 

all 

for 

r»  a n 

then 


[5,  P.  33]  and  [35]. 

U a J f - cV,.  If  all  D f arc  7-Lipschitz  continuous,  then 

k * y **  0 v 

are  7-Lipschitz  continuous;  furthermore,  v U - f 

f 7-Lipschitz  continuous  [43 ,^§30] • From  those  results,  it 
be  seen  that  ^ u =>  f where  U =>  | if  - dv*  . It  0 - En. 
it  is  necessary  to  assume  f = 0(r~2)  nr  a similar  condition. 
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